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Abstract 

We present a search for the decay B — * vr 7r° by the BABAR experiment at the PEP-II asymmetric- 
energy B-factory at SLAC. Using approximately 88 million BB pairs collected between 1999 and 
2002, we place a 90% confidence level upper limit on the branching fraction of 

B(B° -» ttV) < 3.6 x 10~ 6 . 

This result is preliminary. 
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1 Introduction 



The study of B meson decays into charmless hadronic final states plays an important role in the 
understanding of CP violation in the B system. Measurements of the CP- violating asymmetry in the 
B° — > 7r + 7r~ decay mode can provide information on the angle a of the Unitarity Triangle. However, 
in contrast to the theoretically clean determination of the angle f3 in B decays to charmonium final 
states 0, §], the extraction of a in B° — * ir + ir~ is complicated by the interference of tree and penguin 
amplitudes with different weak phases. The time dependent CP-violating asymmetry in B° — > 
7r + 7r~ is proportional to sin2a e ff. Assuming an isospin relation ||, |a e ff — a\ may be determined 
from the branching fractions B(B ± t^tt ), B(B° vr+vr-), B{B° vr+vr-), B(B° -> 7r°vr ), 
and B(B° — > 7r°7r°). Alternatively, a bound on a e g- — a may be found from the ratio B(B° — > 
7r 7r°)/^(5 =l= — > vr^vr ), using the average of B° and 5° branching fractions In this paper, we 
report on a search for the decay B° —* 7r°7r°. Here and throughout this paper B° —* 7r°7r° is meant 
to include both B° and B° decays. 

2 The BABAR Detector and Dataset 

BABAR is a solenoidal detector optimized for the asymmetric beams at PEP-II and is described in 
detail in Ref. ||. Charged particle (track) momentum and direction are measured with a 5-layer 
double-sided silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH) embedded in a 1.5 
T superconducting solenoidal magnet. Neutral cluster position and energy are measured by an 
electromagnetic calorimeter (EMC) consisting of 6580 CsI(Tl) crystals. The photon energy reso- 
lution is a E /E = (2.32/E( GeV) 1 / 4 © 1.85)%, and the angular resolution is a 8 = 3.87°/ v / S( GeV). 
Charged hadrons are identified with a detector of internally reflected Cherenkov light (DIRC) and 
specific ionization in the tracking detectors. The instrumented magnetic flux return (IFR) detects 
neutral hadrons and identifies muons. 

This search uses (87.9 ± 1.0) x 10 6 BB pairs from approximately 81 fb 1 of data at the T(4S) 
resonance (on-resonance) , and approximately 9fb _1 of data at 40MeV below the T(AS) resonance 
(off-resonance), collected with the BABAR detector from 1999 through 2002. The PEP-II collider is 
operated with asymmetric beam energies, corresponding to a boost for the T(4S) of /?7 = 0.55. 

3 Event Selection 

BB events are selected using track and neutral cluster content and event topology. Events are 
required to have either three or more well measured tracks from the interaction point with transverse 
momentum pr > 0.1 GeV/c and polar angle in the lab frame 0.41 < #i a b < 2.54 rad, or two or 
fewer such tracks combined with two or more neutral clusters with center-of-mass (CM) energy 
Ecm > 0.5 GeV and one or more additional neutral clusters with laboratory energy > 0.1 GeV. 
Backgrounds from lepton pair events are removed by requiring that the ratio of the second to zeroth 
Fox- Wolfram moment be less than 0.95 and the event sphericity be greater than 0.01. The principal 
background is from the e + e~ — ► qq process (q = u, d, s, c), when both quark jets contain a 7r° which 
combine to mimic a B decay. This background is suppressed by requiring that the cosine of the 
angle between the sphericity axis of the B candidate and the sphericity axis of the remaining tracks 
and neutral clusters in the event satisfy | cos 9s\ < 0.7. 
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4 Candidate Selection 



Candidate it mesons are formed from two neutral clusters with E > 0.03 GeV whose transverse 
energy profile in the EMC is consistent with that of a single photon. The centroid of the two 
clusters must be separated by at least one EMC crystal. To reduce the background from false 7r° 
candidates, the cosine of the angle between the 7 momentum vector in the ir° rest frame and the 
7T° momentum vector in the lab frame is required to satisfy | cos# 7 | < 0.95. The invariant mass of 
the two photons is required to be within ±3cr of the 7r° mass. 

B — > 7r°7r° candidates are formed from pairs of ir° candidates. The remaining background is 
from qq events that have a spherical topology and pass the | cos#s| requirement, and — > p^ir 
decays in which the tt^ is emitted nearly at rest in the B frame. No other B decay produces a 
significant background for B — > The B^ — > p^ir decay mode has not been observed; the 

limit on its branching fraction is 13(3^ — > p^tr ) < 4.3 x 10~ 5 at 90% CL @. Both backgrounds 
are separated from signal by using the kinematic constraints of B mesons produced at the T(4S). 

The first kinematic parameter is a beam-energy substituted mass tues = \J E% — Vb > where Eb = 
(s/2 + pi ■ ps)/Ei; yfs and Ei are the total energy of the e + e~ system in the CM and laboratory 
frames, respectively, and pi and pb are the momentum vectors in the lab frame of the e + e~ system 
and the B candidate, respectively. The second kinematic parameter is AE = Eb — \/s/2, where 
Eb is the B candidate energy calculated in the CM frame. In B° — > 7r°7r° events the m^s and 
AE resolution are predicted to be approximately 3.8MeV/c 2 and 80MeV, respectively, based on 
simulation. 

There are 3020 candidates with m^s > 5.2 GeV/ c 2 and \AE\ < 0.2 GeV which are used in this 
search. The B° — > 7r°7r° signal efficiency is evaluated with a GEANT4 based detector simulation 0. 
The efficiency to separate closely spaced photons in the EMC is measured using — > Tr^Tt ^ 
and decays, and uncertainty in this efficiency dominates the error in the signal 

efficiency. The B° — > ir°ir° efficiency is (16.5 ± 1.7)%. 

The B^ — > p^ir background is reduced by removing candidates in which the omitted is 
identified. Tracks that are not identified as leptons or kaons, and that are not part of a reconstructed 
Kg — > 7r + 7r~, A — > pir, or 7 — > e + e~ candidate, are used. The track that has a tt^tt invariant mass 
and mES of the tt^tt ?: combination most consistent with the p mass and B^ — > p^Tt hypothesis 
is selected. A cut is applied on a linear combination of the vr^vr invariant mass and the AE 
of the 7r ± 7r 7r° combination which removes roughly 50% of B^ — > p 1 * 1 ^ , with 93% efficiency for 
B° -> 7r°7T . Only (0.40 ± 0.04)% of B ± p ± ir° decays remain after all cuts. 

The qq background that remains after all cuts is further distinguished from signal using a 
Fisher discriminant Tt that combines energy flow and B flavor tagging variables. The energy flow 
variables are Lq = Y^iPii an d L2 = YliPi x \ (3 cos 2 (6^) — 1), where the sum is over all tracks and 
neutral clusters in the event except the daughters of the B — > 7r°vr candidate. Here 9i is the angle 
with respect to the thrust axis of the B candidate and pi is the momentum magnitude, both in 
the CM frame. The B flavor tagging variable is a quality index which classifies the lepton, charged 
kaon, and slow pion vr^ ow (from the decay D* ± — > D°tt^ ow ) content of the event. The quality index 
is ordered by the degree of background rejection. The leptons, charged kaons, and slow pions are 
selected and the events are classified with the B flavor tagging algorithm described in Ref . |Q] . The 
coefficients of J~t are optimized using Monte Carlo simulation of signal and qq background. 
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Figure 1: The distribution of the Fisher discriminant Tt from a fully reconstructed B — > D^nir 
data sample (open circles), and from off-resonance data and on-resonance ?t?es sidebands (filled 
squares). The triple Gaussian parameterizations used in the likelihood fit for B — > 7T TT Si gnal 
(dotted line) and qq background (solid line) are also shown. 



5 Unbinned Maximum Likelihood Fit 

The number of B° — > 7r 7r° events is determined by an unbinned extended maximum likelihood fit 
to 77iES) AE 1 , and Tt- The probability V% a-i) for a given hypothesis is the product of probability 
density functions (PDFs) for each of the variables Xj = (?t?es, AE,Tt) given the set of parameters 
dij. The likelihood function is given by a product over all events N and three signal and background 
components: 

' 3 

.i=i 

The iii are the number of events in each of the three components: B° — > 7r°7r° (n^o n o), B^ — * p^it 
(n p7T o), and qq (n q q). Monte Carlo simulations are used to verify that the fit is unbiased. 
The tties PDF for qq is parameterized by a threshold function || 

/(m E s) = mEsV 1 ~ {m ES /m ) 2 ex.p {-^(1 - (m E s/"io) 2 )}, 

where mo is the average CM beam energy, and £ is found from a fit to on-resonance data with 
| cos 9s | > 0.9. The AE PDF for qq is parameterized by a quadratic function with coefficients found 
from a fit to both on-resonance data in the 7ties sidebands and off-resonance data. The Tt PDF 
for qq is the sum of three Gaussians and is also found using both rrtES sideband and off-resonance 
data, as shown in Fig. |l]. The tbes and AE PDFs for signal and B^ — ► p^n background are 
found from Monte Carlo simulation. The B — * and B^ — > /j ± 7t° tties and AE 1 variables are 
correlated, so a two dimensional PDF derived from a smoothed Monte Carlo distribution is used. 
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The Tt PDFs, shown in Fig. |], for both B° — > 7r°7r° and — ► /r^vr are parameterized as the 
sum of three Gaussians and are found from a sample of fully reconstructed B° — > D^nir events, 
with n = 1, 2, or 3. 

The result of the fit is n w o n o = 23lg° and n q q = 2990 ± 55 events. These statistical errors 
correspond to the point at which log£ changes by one half. The number of B^ — » p^ir events is 
fixed in the fit to n p7T o = 8.4, based on the central value from Ref. f|] of B(B ± —> p ± 7r°) = 2.4 x 1CT 5 
and our estimated efficiency. The distributions of ?7?es> AE, and Tt are shown in Fig. |2| after a 
cut on the probability ratio 

~2 _ n 7T o n o'P 7r o 7r o 

n 7r o 7r o'P 7r o 7r o + n p7T oV p7T o + n q qVqq 

Here the V{ are products of the PDFs for the two other variables, and the n« are the central values 
from the fit. The cut is optimized by maximizing the ratio 



y/n^o^e^o + n p7r oe p7r o + n q -q(iq-q' 

where is the efficiency of the cut. The efficiencies for the ttt-es distribution are 20%, 12%, and 
0.8% for the B° —* 7r°7r°, B± — > /O^vr , and qq components, respectively. The PDF projections for 
each of the fit components, scaled by the appropriate e^, are also shown in Fig. ||. 

The results from the likelihood fit are compared to an analysis that simply uses the number of 
candidates satisfying the requirements 5.260 < tties < 5.289 GeV/c 2 , —0.16 < AE < 0.10 GeV, and 
J~t < —0.20. These cuts were chosen in advance by maximizing the ratio 

N TT<->7rO 

^N 7T 0^+N p7T0 +Nqq , 

where N is the number of events from each source that satisfy the cuts. There are 89 events 
satisfying these requirements. The number of background qq events was determined by scaling the 
number of events with 5.20 < m ES < 5.26 GeV/ c 2 , -0.16 < AE < 0.10 GeV, and T T < -0.20 by 
the appropriate factor given the threshold function describing the w-es distribution. The number 
of background B^ — > p^ir events was estimated using the efficiency from the simulation. We 
find N^o^o = —6 ±11 (stat.). Using simulations based on our PDFs, and assuming a flat prior 
distribution for B(B° — > 7r°7r°), we estimate that there is a 2.5% probability to observe 89 or fewer 
events given the central value of our likelihood fit. 



6 Systematic Uncertainties 

We have estimated the systematic uncertainty in the likelihood fit by varying the PDF parameters 
by their statistical errors, by using different parametrizations, and by varying the B^ —> p^ir 
branching fraction. In each case the likelihood fit to the data is repeated and the change in n^o^o 
is used as the systematic uncertainty. The systematic errors are listed in Table |[ The dominant 
systematic uncertainty is due to the statistically limited sample of data used to parameterize the Tt 
PDF for qq. Since the parameters in the triple Gaussian are highly correlated we transform to the 
uncorrelated parameter space and vary the uncorrelated parameters by ±lcr. The fit is repeated for 
each la variation of the uncorrelated parameters, and the positive and negative changes in n^o^o 
are separately summed in quadrature. The fit is also repeated using an interpolated histogram 
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Figure 2: Results from the maximum likelihood fit. The distributions for a) m,Es> b) AE, c) Tt 
are shown, for candidates satisfying an optimized cut on the probability ratio 1Z. Also shown 
are the PDF projections for signal (dotted line), — > /9 ± 7r° (dot-dashed line), qq background 
(dashed line), and the sum (solid line). These plots do not represent the full information used in 
the maximum likelihood fit, but only a subset of the data. The ratio — log {C/ C max ) is shown in 
d) (solid line) and with statistical errors only (dashed line). 



as the qq Tt PDF, with a change of A(n 7r o 7r o) = —1.1 events. The fit is repeated for values for 
the qq m-Es shape parameter of £ = 24.3 ± 1.3, based on the change in £ as a function of cos 9s- 
The qq AE parameters are varied by their statistical errors. The EMC energy scale is varied by 
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±10.4 MeV based on the statistical error in the mean of AE in the B ± -> /i ± 7r° analysis g, and the 
B° — > 7r°7r° A £7 PDF is changed accordingly. The — ► p^vr veto cut is varied and the changes 
taken as a systematic error. Lastly, the B^ — > p 1 ^ branching fraction is varied from 1.2 x 10~ 5 
to 4.3 x 10" 5 . 

Table 1: Systematic errors on the number of B — * vr vr events in the maximum likelihood fit. 
A±(n 7T o 7T o) are the positive and negative change in the number of signal events from the likelihood 
fit for each systematic source. 



Systematic 


A + (n 7r o 7r o) (events) 


A_(n 7r o 7r o) (events) 


qq J~t PDF parameters 


+7.5 


-2.4 


qq Tt PDF functional form 


+1.1 


-1.1 


qq m E s PDF 


+1.2 


-1.1 


qq AE PDF 


+1.0 


-0.2 


B° tt ^ AE 


+0.8 


-1.1 


B — > p^vr cut variation 


+1.3 


-1.3 


S ± — ► p^vr branching fraction 


+1.6 


-1.9 


Total systematic error on n^o^o 


+8.1 


-3.8 


Efficiency systematics 


10.1% 


-10.1% 


Total systematic 


+8.4 


-4.4 



We calculate the significance of the result, defined as s = — 2 log (£(n^o^o = 0)/C max ), and 
the 90% CL upper limit. The upper limit is evaluated by finding n a where 

C(n)dn = Q g 
f™£(n)dn 

Systematic errors are included in the following way. For the significance, we repeat the fit using the 
changes in qq Tt parameterization, qq /ties parameterization, and B^ — > p ± 7r° branching fraction 
which cause n^o^o to decrease. The — log (£/£ max ) function is shown in Fig. ^d, along with the 
same function before systematic errors are included. The significance of the result is s = 2.5a. The 
systematic errors are included by adding the total systematic A + (n 7r o 7r o), in Table [l|, to nF£ ' . We 
find n^o^o < 46 events at 90% CL. 

7 Results 

To convert the number of events n^o^o into a branching fraction we use 

B (B° -> ^vr ) = n *°"° . 

• N BB - 

N bE = (87.9 + 1.0) x 10 6 is the number of BB pairs in our data sample and the efficiency is e^o^o = 
0.165 + 0.017. The central value of the likelihood fit is B(B° vr 7r°) = (1.6jlg;£(stat.) ±£|(syst.)) x 
10~ 6 . To calculate the branching fraction upper limit we decrease e^o^-o and N BB by one a. The 
upper limit on the branching fraction is 

B(B° -> 7r°7r°) < 3.6 x 10~ 6 at 90% CL. 
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These results are preliminary. The upper limit may be combined with our measurement of the 
branching fraction B(B ± -» v^vr ) = (5.5 ± 1.0 ± 0.6) x 10~ 6 |] to bound the ratio B(B° -> 
Tr°TT°)/B(B ± — ► tt^tt ). Treating the systematic uncertainties in the same way as for the B(B° — > 
7T 7r°) upper limit, and removing correlated systematic uncertainties, we find B(B° — * ir°ir°)/B(B ziz - 
■k^tt ) < 0.61 at 90% CL. Assuming the isospin relations for B — * inr Q this corresponds to an 
upper limit of \a e s — a\ < 51° at 90% CL. 
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